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Cosmologists are exploring two possible sets of explanations for the remarkable observation of cos- 
mic acceleration: dark energy fills space or general relativity fails on cosmological scales. We define 
a null test parameter e{k, a) = n~'^dlnD/dlna— 1, where a is the scale factor, D is the growth rate 
of structure, f2m(a) is the matter density parameter, and 7 is a simple function of redshift. We show 
that it can be expressed entirely in terms of the bias factor, b{a) , (measured from cross-correlations 
with CMB lensing) and the amplitude of redshift space distortions, f3{k,a). Measurements of the 
CMB power spectrum determine QmoHQ. If dark energy within GR is the solution to the cos- 
mic acceleration problem, then the logarithmic growth rate of structure din D/ din a = O^. Thus, 
e(fc, a) = on linear scales to better than 1%. We show that in the class of Modified Gravity models 
known as f{R), the growth rate has a different dependence on scale and redshift. By combining 
measurements of the amplitude of f3 and of the bias, b, redshift surveys will be able to determine the 
logarithmic growth rate as a function of scale and redshift. We estimate the predicted sensitivity 
of the proposed SDSS III (BOSS) survey and the proposed ADEPT mission and find that they will 
test structure growth in General Relativity to the percent level. 
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I. INTRODUCTION 

In General Relativity (GR), there are four variables 
characterizing linear cosmic perturbations: the two gravi- 
tational potentials ^' and <i>, the anisotropic stress, a, and 
the pressure perturbation, 6p. All of these variables can 
depend on the wave number, k, and the expansion fac- 
tor a. We focus initially on models with no dark energy 
clustering or pressure and discuss these effects later in the 
paper. We assume scalar linear perturbations around a 
flat FRW background in the Newtonian gauge. 



(1 + 2'^)dt^ + a^{l + 2$)da;^ 



(1) 



and work in the quasi-static, linear approximation, which 
is valid for sub-horizon modes still in the linear regime. 
The evolution of perturbations is described by the conti- 
nuity, Euler and Poisson equations {i.e. [1]): 



VL + K, 



kuV^ 



aH 



(2) 
(3) 
(4) 



With the assumption of no anisotropic stress, <i> = — ^E", 
these equations can be combined to derive the equation 
of motion for the growth factor D, defined as D{k,z) = 
5yn{k,z)/5yn{k,z = 00): 



D" + (2 - 



H 



D' 



1p 



■praD = 0; 



(5) 



where H = a/a is the Hubble function, and a prime de- 
notes derivative with respect to In a. 
From Eq. (5) one can infer the two key features of GR 
with smooth dark energy. First, the growth factor is 
exactly determined once the Hubble function H{a) is 
known; and second, since none of the coefficients is a 
function of scale, the growth factor is scale- independent. 
Therefore, for a given expansion history, one can test GR 
in two ways: checking that the theoretical solution of Eq. 
(5) agrees with observations, and testing the hypothesis 
of scale-independence [2-11] . The growth rate of struc- 
ture in GR is well approximated by fi^, where the fitting 
function 7(2) ~ 0.557—0.022: is accurate at the 0.3% level 
[12]. We define a function that tests the growth rate of 
structure and can be directly related to observables: 



e(fc, a) 



_Ci-fia)d}nD_ 
" " dlna 



a^''H{a)^'' dlnD 
{nm^oH^y dlna ^ ' 



The combination flm.oHQ can be constrained via CMB 
measurements; it is currently known to within the 5% 
level fro the WMAP5 data [13], with an expected gain of 
a factor ~ 4 from the upcoming satellite CMB mission 
Planck [14]. The solid line in Fig. 1 shows e(a) in 
GR with dark energy: regardless of the details of the 
dark energy model, e(a) ~ in the linear regime. By 
measuring this quantity, we can characterize deviations 
from General Relativity. 



II. F(R) THEORIES IN THE PPF FORMALISM 
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whose action is written as 



S = J d^x^g[R + f{R)) + 5„ 



(7) 



Sm is the action of standard matter fields. It has been 
noted long ago [15, 16] that models where f{R) is an 
inverse power of the Ricci scalar can give rise to late- 
time acceleration; however, many of those models have 
been shown not to be cosmologically viable due to grav- 
itational instability [17]. 

Recently, Song, Hu and Sawicki [18] introduced an ef- 
fective parametrization for f{R) theories which does not 
rely on any particular model and is able to discard models 
suffering from instabilities. The cosmological evolution is 
obtained by fixing the expansion history to match that 
of a dark energy model, for which we assume a constant 
equation of state w: 



-3(l+tu)^ 



(8) 



Such requirement for H{a) translates into a second-order 
equation for f{R), which can be solved numerically. Of 
the two initial conditions of this equation, one can be 
fixed requiring that f{R)/R ^ at early times {i.e. for 
large R) in order to recover GR; the second defines a one- 
parameter family of curves which all generate the given 
H{a). Such parameter is conveniently chosen as (also see 
[19]): 



h 



H 



RR. n/ 

1 + fR H 



(9) 



Here and fuji represent the first and second derivative 
of / with respect to R, respectively. GR is represented by 
the special case i?o = 0, so that Bq effectively quantifies 
the deviation from GR at the present time. Furthermore, 
the gravitational stability condition is easily established 
as Bq > 0. 

The additional degrees of freedom of the f{R) gravity 
introduce modifications in the Poisson equation and in 
the relation between the two gravitational potentials [1] , 
which now read: 



and 



^{g{k,a)-l) — i^-^). 
all 



(10) 



(11) 



The equation for the growth factor in MoG takes the 
form: 

p. , , H' 4^Geff(fc,z)(g(fc,z)-l) 

D +{2 + —)D + — pmD = 0. 

(12) 

Given H{a), a MoG model is not completely defined un- 
less G'off(fc, a), the effective Newton's Constant, and the 
metric ratio g{k, a) = ($ -I- 5')/(<I> — ^), are known, con- 
trary to the GR case. For the class of f{R) theories 
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FIG. 1: The behavior of t{k,a) = Q.'^ dlnD / dlna - 1 in 
GR (solid line) and in /(ii) models, as a function of Bo and 
k. Growth is enhanced for Bo 7^ and at smaller scales in 
alternative theories. In GR, e(a) = 0. 



under study, and in the sub-horizon, linear regime, Hu 
and Sawicki [20] have provided a fit to Gofr and g as 



Geff(fc,a) 
g{k,a) 



Gcfr(a) = 



G 



N 



l + ./fi' 



(13) 



gsH{a) - \{0.nVB{k/aH)f 



1 + {0.71V B{k/aH)) 



where B is the function appearing in Eq. (9), and gsH{o-) 
is the super-horizon metric ratio (see [18, 20] for details). 
The key result is that the f{R) models all predict an 
enhancement in the growth rate of structure. In fact, 
any positive value of Bq gives rise to a negative fa, 
so that the effective Newton constant is larger with re- 
spect to GR. Moreover, both terms in g{k, a) have a neg- 
ative sign, which induce further enhancement of mat- 
ter clustering, as can be seen from Eq. (12). Since 
Gcff(fc, z) — Gn I {1 + Ir) does not depend on k, the scale 
dependence of the growth factor can only arise from the 
second term of Eq. (13). On large scales, the dominant 
term in the metric ratio is (a), while for increasing 
values of k, the second term in the expression for g{k^ a) 
becomes important and tends to the constant value of 
— 1/3 for large k. The scale of the transition from scale- 
free to scale-dependent growth factor is k/aH ~ B'^^^; 
due to the asymptotic behavior oi g{k,a), the growth dif- 
fers significantly from GR even for models with very small 
values of i3o, on sufficiently small scales. A mechanism 
to restore GR on scales of the galaxy and smaller is dis- 
cussed in [20, 21]. Fig. lshowsdln(D)/dln(a)f}„i(a)''-l 
for a few different values of Bq and k/aHo. For the GR 
case, Bq — 0, £{k, a) — 1 ~ with no scale dependence. 



III. MEASURING THE GROWTH OF 
STRUCTURE 

Peculiar velocities displace galaxies along the line of 
sight in redshift space and distort the power spectrum of 
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galaxies observed in redshift space. This effect is known 
as linear redshift space distortion and was first derived 
by Kaiser [22]. In redshift-space, the power spectrum 
is amplified by a factor (1 + over its real-space 

counterpart, 



(14) 



where P^(k) and P{k) are redshift and real space power 
spectra, respectively, /Xk = 5 • fc is the cosine of the angle 
between the wavevector k and the line of sight z, and (3 
is the linear redshift-space distortion parameter defined 

as 



I dlnD 
b dlna 



(15) 



b is the linear bias, which we assume to be independent 
of scale. 

Galaxy redshift surveys can be used to directly mea- 
sure (3, and, if the bias is known, the growth rate of per- 
turbations. In fact, the redshift space power spectrum 
can be decomposed into harmonics, whose relative am- 
plitude depend on the growth rate of structure through 
/3. The possibility of using such dependence in order to 
constrain dark energy properties has been explored in 
[23]; here we focus on the measurements of the scale de- 
pendence of /3 as a smoking gun of Modified Gravity. 
We assume that /3 is obtained through the ratio of 
quadrupole to monopole moments of the redshift power 
spectrum [24] 



P2(fc) _ §/3 + 



Po(fc) l + |/3 + i/32' 



(16) 



and use the prescription in [25] to get the errors in the 
above quantities: 



1/2 



Vk [/d3rn2(r)V'2(r)]' 



(17) 

where n(r) is the mean galaxy density, ip is the weight 
function, Vk is the volume of the shell in A;-space, and 
the index "i" assumes the values and 2 for monopole 
and quadrupole, respectively. 

The linear bias for a population of large-scale structure 
tracers can be estimated by cross-correlating the line-of- 
sight projected density of the tracer with a convergence 
map reconstructed by CMB Icnsing techniques, and com- 
paring the resulting signal with theory. The weak lensing 
potential responsible for lensing the CMB can be written 
as the line-of-sight integral [26], 



where ^^(77) is the comoving angular diameter distance 
corresponding to the comoving distance rj, and rjo is 



the comoving distance to the last scattering surface. A 
quadratic combination of the measured CMB tempera- 
ture and polarization [27 29] provides an estimator of 
the convergence field, k = ^V^*/). In this study, we have 
used the prescription of [27] to compute the expected 
noise power spectrum, N^^'^ , corresponding to the recon- 
structed convergence field by cross-correlating k with the 
projected fractional overdensity of the tracer, 

E(n)= j dvW{v)b6m{vn,v) (19) 

where W is the normalized tracer distribution function 
in comoving distance. We measure the cross-correlation 
spectrum: 



c: 



-b CtmH, 



■2 ^ 



{rf) dA' dA{r])dA{rio) 



(20) 

where P{k,ri) is the matter power spectrum at the co- 
moving distance r] and we have related the wavenumber 
k to the multipole i via the Limber approximation [30]. 
The signal-to-noise ratio for such a cross-correlation can 
be estimated as [31], 



-) =/.k,i:(2^+i) 



(21) 

where /sky is the fraction of sky over which the cross- 
correlation is performed. For tracer counts the noise is 
Poisson, and the power spectrum is given by, Nf''^ = 1/h 
where h is the number of tracer objects per steradian. 

Since the signal is proportional to the bias, b, the ex- 
pected error on b can be written as 



M/bc^l/{S/N). 



(22) 



We consider throe present and forthcoming redshift 
galaxy surveys: the SDSS LRG sample [32] , its extension 
BOSS-LRG [33], which we divide in two redshift bins, la- 
beled as BOSSl and B0SS2, and the proposed survey 
ADEPT [34]. Specifics of each experiment are listed in 
Table I. In all cases we assume that /3 and b do not change 
significantly with redshift within a survey, so that the ob- 
served quantity is (}(k, Zc), where Zc is roughly the central 
redshift of the survey. As a direct comparison the capa- 
bilities of the three galaxy surveys under examination, we 
show the real space matter power spectrum, normalized 
to the SDSS LRGs median redshift, with its crrorbars 
in Fig. 2. We also consider three possible CMB experi- 
ments: a PLANCK-hke CMB experiment with 65% sky 
coverage and temperature and polarization sensitivities 
of 28 /iK-arcmin and 57 /xK-arcmin, respectively; a next 
generation CMB survey based on using a camera simi- 
lar to that on ACT or SPT with a polarimeter and a 
~ 3 years observing program (labeled PACT) with 65% 
sky coverage and temperature and polarization sensitiv- 
ities of 13 /iK-arcmin and 18 /iK-arcmin and an ideal 
polarization experiment (labeled IDEAL), with 65% sky 
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FIG. 2: Errors on P{k), normalized to the SDSS-LRG median 
redshift (z=0.31) for all surveys. 
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TABLE L Predictions for the errors on bias from the cross- 
correlation studies described in the text. For each combina- 
tion of experiments, we display the number of galaxies per 
square degree (n); the area of overlap (A), the signal-to- noise 
with which the cross correlation of tracer surface density with 
CMB-lensing can be extracted, (S/N), and the percentage er- 
ror in the bias, b, for the tracer. 



coverage and temperature and polarization sensitivities 
of 1 /.tK-arcmin and 1.4 /iK-arcmin, respectively. The 
expected results from cross-correlation with the ADEPT 
and BOSS surveys, and the SDSS LRG are displayed in 
Table I. 



IV. RESULTS 

We show our main results in Fig. 3. Errorbars are com- 
puted using Planck as the complementary CMB lensing 
survey for SDSS LRG and BOSS, and PACT for ADEPT. 
We summarize the current status of /3 and bias measure- 
ments in Table II, and add the expected errorbars on e 
at two different scales, k/h = 0.05 Mpc~^ and k/h = 0.2 
Mpc~^, from our analysis, for comparison. 



The errors on bias are scale-independent and vary from 
17% for SDSS LRG to ~ 1 % for ADEPT. Errors on (3 
depend on scale; we bin our simulated data in bins in k 
space of width Alnfc = 1/e, and see that for all surveys 
errorbars decrease as we go to smaller scales. The corre- 
sponding constraints on e are as strong as a few percent 
for BOSS, and of the order of 1% for ADEPT. 

The combined effect of smaller errors and of the asymp- 
totic behavior of the growth factor, which induces a large 
deviation of e{k,a) from its GR value on small scales, 
is that redshift galaxy surveys are more sensitive to the 
small-scale modification of gravity than to the large-scale 
one. Ultimately, the smallest observable value of Bq will 
not be set by the capabilities of the survey, but by the 
breakdown of the linear regime assumption. Assuming 
k ~ 0.05 Mpc~^ as an upper limit, and for values of aH 
corresponding to redshifts between and 2, the smallest 
Bq inducing scale-dependent growth is Bq ~ lO"**. Such 
value is within reach of ADEPT; future experiments that 
detect redshifted 21 centimeter emission could probe even 
larger values of k/aH in the linear regime. 



V. CONCLUSIONS AND DISCUSSION 

We have built a null test parameter for General 
Relativity, e(fc,a), based on the consistency between 
expansion history and structure growth expected in 
GR. Such parameter can be expressed in terms of the 
combination f2om^^, probed by the CMB experiments, 
the linear matter perturbations growth factor, probed 
by redshift galaxy surveys, and the linear bias, probed 
by cross-correlation of the two. 

We have predicted the achievable precision in the 
measurement of e(fc, a) for three redshift galaxy surveys, 
SDSS LRG, BOSS and ADEPT, together with Planck 
and a possible future CMB experiment, PACT. We have 
interpreted such result in the context of a one-parameter 
family of modified gravity theories, known as f{R), 
which can give rise to cosmic acceleration. In such 
models, the matter clustering is enhanced on all scales 
with respect to the GR case, and the enhancement 
is largest on small scales. We concluded that the 
peculiar signatures of the f{R) theories will be definitely 
detectable with a survey like ADEPT. 
More generally, any detection of deviation of e from 
zero that was not due to some observational system- 
atic would be a signature of truly novel physics with 
enhanced growth, pointing either to non-GR physics or 
to unexpected properties of dark energy: dark energy 
models with a non-zero sound speed are characterized 
by an oscillatory behavior of the growth [35] , and scalar 
field dark energy suppresses growth on large scales [36]. 
Similarly, massive neutrinos suppress e(fc, a) on scales be- 
low the neutrino free streaming scale (see [37] for review). 

We warmly thank E. Aubourg, C. Hirata, W. Hu, R. 
H. Lupton, M. A. Strauss and L. Verde for useful sugges- 
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FIG. 3; e(fc, z) for the four surveys, as a factor of Bo and k. Total errorbars around the ACDM case are shown in black; the 
smaller red errorbars are from bias only. 
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